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We present a numerical integration of the cosmological scalar perturbation equations in warm 
inflation. The initial conditions are provided by a discussion of the thermal fluctuations of an inflaton 
field and thermal radiation using a combination of thermal field theory and thermodynamics. The 
perturbation equations include the effects of a damping coefficient T and a thermodynamic potential 
V. We give an analytic expression for the spectral index of scalar fluctuations in terms of a new 
slow-roll parameter constructed from T. A series of toy models, inspired by spontaneous symmetry 
breaking and a known form of the damping coefficient, lead to a spectrum with n s > 1 on large 
scales and n s < I on small scales. 



I. INTRODUCTION 



Inflation is the most successful idea which we have available for explaining the large scale structure of the universe. 
Inflationary dynamics can be realized in two distinct ways. In the original picture, termed supercooled inflation, the 
universe rapidly supercools during an inflation phase and subsequently a reheating period is invoked to end inflation 
and fill the universe with radiation UJ, |3i S Ij] • I n the other picture, termed warm inflation, dissipative effects are 
important during the inflation period, so that radiation production occurs concurrently with inflationary expansion. 

The idea of warm inflation was influenced by a calculation of the friction term in the inflaton equation of motion 
by Hosoya and Sakagami The magnitude of the damping term suggested the possibility that it could be the 
dominant effect prolonging inflation G, 7, 8. 9]. The significance of dissipation to inflation was independently realized 
in [T(il |. but this time with a clearer picture of the associated inflation dynamics, and named warm inflation. This 
and associated works [TTL Il2| also outlined the nonequilibrium thermodynamical problem underlying this picture, 
which has subsquently has been studied in greater detail |13L fl4L Il5t Ha . H71 Ha . 119] . These works confirmed that at 
high temperature the damping was proportional to the relaxation time of the radiation. In addition to these studies 
of warm inflation dynamics, several phenomenological warm inflation models have been discussed in the literature 

The density fluctuations in warm inflation arise from thermal, rather than vacuum, fluctuations [|| l25l l26j | . These 
have their origin in the hot radiation and influence the inflaton through a noise term in the equations of motion 

mm. 

The development of cosmological perturbations in warm inflation has been investigated by several authors. Taylor 
and Berera |27| have analysed the perturbation spectra by matching the thermally produced fluctuations to gauge 
invariant parameters when the fluctuations cross the horizon. This technique should be good for order of magnitude 
estimates. For more accurate treatment of the density perturbations it is nessesary to use the cosmological perturbation 
equations, which can be found in the literature psl l29l l30l|. 

The main focus of the present work is to solve the perturbation equations numerically. We will also set up the equa- 
tions to take into account some effects which have previously been ignored. We introduce a new slow-roll parameter 
and include modifications to the equations when the damping term and the potential depend on temperature. 

Other improvements have been made in the treatment of thermal fluctuations. The influence of the cosmological 
expansion on the inflaton fluctuations is considered in detail. This gives a clearer and more accurate picture of the 
development of the thermal fluctuations than previously. Thermal fluctuations in the radiation, which lead to entropy 
perturbations, are also analysed. Entropy perturbations are always present during warm inflation and react back on 
the curvature fluctuations. In the basic model, the entropy fluctuations disappear before inflation ends. 
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Section 2 contains an updated outline of the theory of warm inflation. In section 3 we obtain expressions for the 
thermal fluctuations of the inflaton and the radiation. Cosmological perturbations are discussed in section 4. The 
results of numerical solutions for the cosmological perturbations using the thermal fluctuations as initial conditions 
are compared with observations in section 5. A summary of the main results and some general comments appear in 
section 6. 



II. WARM INFLATION 



In this section we shall consider the case of a homogeneous inflaton field interacting with thermalised radiation. 
If the inflaton is evolving very slowly then the radiation creates a thermal correction to the inflaton potential |3ll | 
and a damping force [f| . The system of equations describing the system include an equation for the inflaton and an 
equation describing how the energy lost by the inflaton through the damping force is transfered to the radiation. 

The equation describing the effects of radiation damping on the evolution of the inflaton field was first obtained by 
Hosoya || using a form of thermodynamic transport theory. This is applicable as long as the inflaton evolves slowly 
and the radiation remains close to thermal equilibrium. More recent treatments are based on the closed time path 
formalism, and an account of the effects of radiation damping on an evolving field in flat space can be found, for 
example, in a paper by Gleiser and Ramos l.32i. An interesting variation occurs when the inflaton interacts with the 
radiation via an intermediate particle decay |19| . 

When the early universe is in homogeneous expansion with expansion rate H(t), the evolution equation for the 
inflaton field is given by 

<£ + (3# + r)^ + ^ = o, (i) 

where r(</>, T) is the damping term and V((j), T) is the thermodynamic potential. The damping term has a generic 
form given approximately by T ~ g 4 (f> 2 T, where g is a coupling constant. For Hosoya's damping term, r = r(<p,T) 
is related to the relaxation time of the radiation and for the models with an intermediate particle decay r = t(0) is 
related to the lifetime of the intermediate particle. The thermodynamic potential can be expressed in a standardised 
way by [3^| 

T) = -^9*T 4 + \{6m T fp + V (ct>), (2) 

where g*(T) is the effective particle number and (5?71t(</>, T) represents thermal corrections. 

The relative strength of the thermal damping compared to the expansion damping can be described by a parameter 

r, 

Warm inflation is defined to be inflation with large values of r, as opposed to supercooled inflation in which r can be 
neglected. 

The dissipation of the inflaton's motion is associated with the production of entropy. The entropy density of the 
radiation s(</>, T) is defined by a thermodynamic relation in terms of the thermodynamic potential, 

s = -V T - (4) 

The rate of entropy production can be deduced from the conservation of energy-momentum. The total density p and 
pressure p are given by 



P 



U 2 + V + Ts (5) 



2 

1 i2 



P = iP - V. (6) 
Energy-momentum conservation, 

P + 3H(p + p)=0, (7) 
now implies entropy production. Making use of Eq. Q we get 



T{s + 3Hs) = I> 2 . 
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The zero curvature Friedman equation completes the set of differential equations for 4>, T and the scale factor a, 

3H 2 = 87tG(±</> 2 + V + Ts). (9) 

The entropy production has been described in a slightly different way in some previous work on warm inflation 
|lCt [2(| . We can recover an alternative equation in the case when the temperature corrections to the potential are 
negligeable. If we set <5m*r = in Eq. (J2J), then the radiation density p r = AsT/3, and Eq. (jHJ becomes 

p r + 4Hp r = r<p 2 . (10) 

This equation is only valid when Suit = 0. 

The slow-roll approximation consists of neglecting terms in the preceeding equations with the highest order in time 
derivatives 

* = -mrr) <"> 

Ts = r<j) 2 (12) 
3H 2 = 8ttGV. (13) 

Note that Ts is the same order as two time derivatives. Slow-roll automaticaly implies inflation, a > 0. 

The consistency of the slow-roll approximation is governed by a set of slow-roll parameters. Warm inflation has 
extra slow-roll parameters in addition to the usual set for supercooled inflation due to the presence of the damping 
term T. There are four 'leading order' slow-roll parameters, 

-M^y '-ds(W 

The first two parameters are the standard ones introduced for supercooled inflation [34l l35l l36l l37j . Two new 
parameters are required for the <p dependence of the damping term and the temperature dependence of the potential. 
The slow- roll approximation is valid when all of the slow-roll parameters are smaller than 1 + r. In supercooled 
inflation, the condition is tighter and the slow roll parameters have to be smaller than 1. 

As an example, consider the case = V.<j>T = 0. In the warm inflationary regime r > 1, the rate of change of 
various slowly varying parameters are given by differentiating equations I|llll3fl . 

1 dlnH 1 , . 

(15) 



H dt r 

- -;<"-«■ < M > 

1 d In Ts 1 . n % , 

H at r 

In each case, the right hand size of these equations gives the relative sizes of terms neglected in the slow-roll approx- 
imation. Slow-roll therefore requires e <C r, t/ <C r and (3 <C r. 

Inflation ends at the time when a = H + H 2 falls to 0. From Eq. I|15(l . the end of inflation corresponds to the 
condition e = r. Furthermore, the slow-roll equations allow us to rewrite this condition as p r — V, in other words 
inflation ends when the radiation energy density reaches the same value as the vacuum energy. 

Similar arguments apply when T t and V^t are non-zero, for example, 

1 dm</> 1 f A - c 4 c 3c \ 

:V-T—J+T—J-T—J) (18) 



H dt r \ 4 + c 

Note that c = TT^/T is not required to be small and therefore does not enter the list of slow-roll parameters. The 
slow-roll approximation is valid for e r, r/ r, /3<^r and i5<Cr. 



III. THERMAL FLUCTUATIONS 



In this section we shall consider the behaviour of thermal fluctuations during inflation and their influence on density 
perturbations on length scales smaller than the size of the horizon. Thermal fluctuations, if they are present, form 
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the predominant source of density perturbations in warm inflation. We will assume throughout that the radiation is 
close to thermal equilibrium. Surprisingly, this is not a necessary condition for warm inflation but it is an important 
regime. Conditions for thermalisation are model dependent, and are discussed in refs [l3l ITil liH ITflj . 

The evolution of a comoving mode of the inflaton fluctuations during a warm inflationary era can be divided into 
three regimes, depending on the relative effect of different physical processes: 

1. Thermal noise 

2. Expansion 

3. Curvature fluctuations, 

The transitions between these regimes are called freezeout and horizon crossing. Where previous calculations have 
treated each of these regimes seperately, we shall consider the combined effects of the thermal noise and the expansion, 
and then the combined effects of expansion and curvature fluctuations. The main results of this section are equations 
1)370. (|38|l and (|51|l . which provide the initial conditions for the classical evolution equations described in the next 
section. 



A. Inflaton fluctuations 



The behaviour of a scalar field interacting with radiation can be analysed using the Schwinger-Keldysh approach to 
non-equilibrium field theor y |38 l. l39j . According to this approach, a simple picture emerges in which the field evolves 
by a Langevin equation |32l |40| , 

- V 2 cf>{x, t) + T<j>(x, t) + = £(x, t), (19) 

where £ is a stochastic noise source. Correlation functions can be evaluated by probability averages. 

The Schwinger-Keldeysh approach has been studied for curved, as well as flat, spacetime backgrounds 41] . Although 
the particular form of the Langevin equation has not yet been fully investigated, we can appeal to the equivalence 
principle to deduce the form of the Langevin equation on scales smaller than the horizon 26]. 

We shall analyse the inflaton fluctuations on an expanding background over timescales which are longer than the 
microphysical processes, but short compared to the variation of the expansion rate. We also take modes with physical 
length scales larger than the microphysical processes. The comoving wave number will be denoted by k and the 
physical wave number q = ka^ 1 . 

The Langevin equation for the fourier transform of the inflaton fluctuations r\k takes the form |26| ] 

^. + 3ffr) fc + r7) fe + fc 2 a- V (20) 

where ^ is a stochastic noise source. Over the timescales of interest, H and T are constant and a = exp(Ht). As in 
the previous section, we set r = T/3H. 

The correlation function of the noise can be found by converting flat space results |32| to the comoving frame. In 
the high temperature limit T — > oo, the noise source is Markovian, 

(*0> = 2rra- 3 (2^) 3 <5 3 (fc - k')S(t - t'). (21) 

The scale factor appears due to the use of comoving coordinates. (The temperature should be large in comparison with 
the mass of the particles which make up the radiation for this approximation. We shall make the same approximation 
for the damping term later.) 

For the eventual comparison with observational data, we define the power spectra V(k) by 

P^ik^irk^fS^k-k') = (t)) (22) 

V^{k){2nk- 1 ) 3 S 3 {k-k') = (fik{t)ri-k'(t)) (23) 
P^(fc)(2^fc- 1 ) 3 ( 5 3 (fc-fc') - (flk{t)V-k'(t)) (24) 

When comparing to classical perturbation theory it is also useful to define root-mean-square fluctuation amplitudes 

8<f>(k) = |7^(fc)| 1/2 (25) 
ty(fc) = |7V(fc)| 1/2 (26) 



Defined in this way, 5<f> has the same dimensions as cj). 
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The power spectra can be found by solving the stochastic equation (|20|l by green function methods, 



m{t) = ~Jpj z (27) 
where = k/a(t)H and the order of the Bessel functions is 

v=l(l + r). (28) 
It is now possible to use the relation 121|) to obtain the power spectrum, 

- - ! 2rr \j v {z)y v {z') - ^(z'^MlV'V 2 - 2 ^'. (29) 

Similarly, for the time derivatives, 

TWfc) = - (I)' H 2 z 2 2TT £ {J v ^{z)Y u {z') - J v {z')Y v ^{z)} 2 z 2 » ' z a ^ ' dz> . (30) 

The formula for the power spectrum can be simplified at late or early times. For early times, when the physical 
size of the mode a/k is much smaller than the horizon, we can use the z ^> v approximation for the Bessel functions 
to obtain 

V^(k) ~ ka~ l T (31) 
V^{k) ~ k 3 a- 3 T (32) 

Allowing for the transformation from comoving to physical wavenumbers, these results are identical to the thermal 
fluctuations of a field in flat space. The rms fluctuations are proportional to T 1 / 2 , which is typical of the fluctuations 
of a classical object in a radiation field [42| . 

A different approximation holds for later times, when H < ka^ 1 < (TH) 1 ^ 2 , and we can use the z <C v approxima- 
tion for the Bessel functions. The power spectrum simplifies to 

„2 poo 

/~\2 / j ( l\2 J1-2v 



~ —YTz AV Y v {zY \ J v (z'Yz' A - il > + O(z') (33) 
2 Jo 

The root-mean-square fluctuation is approaching the dominant solution z v Y v {z) of the homogeneous equation for rjk 
in this regime as if the noise had no effect. With the warm inflation condition T H, the results further reduce to 



vMk) ~ G) 1/2 f^™ 1/2T (35) 

The cross-correlation reduces to 

V^{k) ~ -(^(fc)) 172 ^™^)) 1 / 2 , (36) 

which also indicates that the effects of the noise term are fading away. 

The transition between the two regimes has been called freezeout \2(\ . and occurs at a particular time tp-(fc) when 
ka" 1 = (TH) 1 / 2 . In warm inflation, the freezeout time will always preceed the hubble crossing time, at which 
fca -1 = H . The evolution of the inflaton becomes increasingly deterministic for times t > tp. 

The fluctuations approach a constant value, 

S<Kk) - {YHf/^l 2 (37) 



~ - (?) ^{m 1 ^ 1 ' 2 . (38) 



7T\l/4 k 2 

rv 

Apart from the numerical factor, the value of <5</> agrees with the results of Berera [2(|, who analysed equation 12U|) 
in the regime ka" 1 > {YHf/ 2 . The sign of 8(f> has been chosen for consistency with the cross correlations (|36|) . It is 
also consistent with the deterministic solution. 
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B. Energy fluctuations 

Now we turn to thermal fluctuations in the energy density in flat space. The classical theory of energy fluctuations 
is based on fundamental principles of statistical physics. For the energy fluctuation AE in a fixed volume V, Lifshitz 
and Pitaevstii ji^ gi ye the formula 

((AE) 2 )v = C V T 2 , (39) 

where C v is the specific heat. This result can also be expressed in terms of the energy density p and the entropy 
density s, 

((Ap) 2 ) v = r'r 3 ^. (40) 

We shall investigate the equivalent result in finite temperature quantum field theory before considering the application 
to warm inflation. 

Consider a scalar field \ m thermal equilibrium at temperature T. Fluctuations e in the energy density are defined 

by 

e = H(x)-{H) T , (41) 

where the energy density operator 

^X 2 + i(Vx) 2 + ^mV, (42) 

and the energy density p r = (H)t is the thermal ensemble expectation value of Ti.. 
The correlation function 

(e(x)e(x')) T = (H(x)H(x')) T - (H) 2 T . (43) 

The simplest way to remove any divergences will be to subtract the zero temperature correllations, which in any case 
are not of interest to us. We therefore define the power spectrum of the radiation V r by 

TAq^irq-'fS^q - q') = {e q {t)e^{t')) T - (s q (t)e^ (t')) . (44) 

As before, we also define a root-mean-square fluctuation 

S Pr (k) = iVriq)^ 2 . (45) 

The expectation values for four field operators can be reduced using the cluster decomposition principle |4flj to give 

d 3 q' 



V r (q) = J ^3 q 3 (-q(q - <?') + d t d V + m 2 ) 2 {G>{q, t)G>(q - q' , t') - G>(q, t)G% (q - q\ t% 



(46) 



where the expectation values of products of two fields is given by the thermal green function, 

( Xq (t)X- q >(t')) - (2n) 2 6 3 (q - q')G>(q,t- t'). (47) 
For a bosonic field in thermal equilibrium, 

G^{q,t)=n{u q )e i "<> t + {l+n{w q ))e- iu »\ (48) 
where uj q = (m 2 + q 2 ) 1 ^ 2 and n{uj) — (e"/ T — l) -1 . The power spectrum reduces to the integral 



V M = I Aufuj^, (H 2 - ? ' w 9' - m V) (2n(w,_ fl 0n(w,/) + n(w,_,0 + n(uy)) ■ (49) 



The classical limit corresponds to q -C T. The integral is then dominated by the region q' 3> q. For a massless 
field, the integral reduces to 

V r {q) =M 3 T Pr , (50) 
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where p r is the energy density. This is in agreement with the thermodynamic result (|40|l for pure radiation if we 
identify q 3 with V -1 . 

In the expanding universe, the inflaton and radiation fluctuations are effectively in equilibrium for times t < tp, 
where tp is the freezeout time mentioned earlier. The total energy fluctuations in this regime are given by equation 
(|40|l . After tp, the radiation and inflaton fluctuations are uncorrelated and equation l|50|) can be used for the 
fluctuations in the energy density of the radiation, 

Spr(k) = (^) 1/2 gl /2 k^a-^\ (51) 

where g* is the effective particle number and we have introduced the comoving wave number k = aq. 
The total energy is defined by equation JSJ. The inflaton energy fluctuations for t > tp are therefore 

5p4 = <j>6<t> + V t t5cj), (52) 

where S(f> and S(f> are given above. 



IV. COSMOLOGICAL PERTURBATIONS 



The thermal fluctuations occuring during warm inflation evolve gradually into cosmological perturbations. The 
thermal noise has a diminishing effect on the development of the perturbations until the freezeout time when it becomes 
insignificant. The perturbations can then be described by a gaussian random field which grows deterministically. 

In warm inflation there are cosmological perturbations in the inflaton field, the radiation and the gravitational field. 
We shall only consider the scalar gravitational mode, which means that there is one degree of freedom in the metric 
perturbations. We shall also use a particular gauge, the zero-shear gauge, in which the scalar metric perturbation has 
the form 

ds 2 = -(1 - 2tp)dt 2 + a 2 (l + 2tp)8 ij dx i dx j . (53) 

The inflaton perturbation will be denoted by 8<j), the temperature perturbation ST and the velocity perturbation v. 
For any random perturbation field 3(2;), the power spectrum is defined by an average 

V g (k)(2TTk- 1 ) 3 5 3 (k - k') = {g k g- k >), (54) 

where g k is the fourrier transform of g{x). The amplitude will be normalised by 

g(k) = \V g (k)\ 1/2 - (55) 

Note that g(k) satisfies the same linear equations of motion as g k . We use this normalisation for the perturbations 
and omit explicit reference to k. 

The total energy density and pressure perturbations are then 

Sp = <j>54> + V t 4,54> + <fi 2 (p + TSs (56) 
Sp = ^S4>-V^8(f> + 4) 2 ^ + s8T (57) 

and the energy momentum flux (p + p)v. 

A com plet e set of perturbation equations can be obtained from the Einstein field equations and the scalar field 
equation [2|| y(j • From the Einstein equation we obtain 

<p + H<p + 4nGk- 1 a(p + p)v = (58) 
3Htp+(3H 2 + k 2 a- 2 )tp- AnGSp = (59) 
tp + AHip + (2H + 3H 2 )tp + AnGSp = 0. (60) 

Perturbations of the scalar field equation give 

S4> + (SH + T)8<j> + <f>{8T) + k 2 a- 2 S(j) + SVtf, + A^tp - Tfap - 2V^tp = 0. (61) 
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In general, we can have T = T((f), T) and V = V((f>, T), 

ST = T^54> + r iT ST (62) 

SVf = V^ + V^ T ST (63) 

Ss = -V^tH-VttST. (64) 

The numerical results of the next section are obtained by integrating these equations with the initial conditions set 
by the thermal fluctuations. 

During warm inflation, the background solution rapidly approaches the slow-rolldown approximation lllll3fl . In a 
similar way, the perturbations also have a slow-roll limit, which begins when their length scales become larger than 
the horizon, k < aH. The <5</>, <p and ip terms become insignificant, and perturbations approach a slow-roll large-scale 
limit 

5cj) CH- 1 ^ (65) 

tp ~ 47rGC(l + r)0 2 (66) 
v Cka^H- 1 (67) 

where C is a constant and r = T/3H. 

The cosmological perturbations can also be described in terms of gauge invariant quantities, and these are useful for 
following the development of perturbations after the end of inflation. The curvature perturbation 1Z [44| , is defined 
by 

1Z = tp- k- 1 aHv, (68) 

The entropy perturbation e is defined by 45] 

e = Sp — c^Sp, (69) 

where c 2 = pj p. In the slow-roll large-scale limit, the curvature perturbation is constant, 1Z ~ C and the entropy 
perturbations vanish. 

The curvature and entropy perturbations are coupled by a second order equation [30| 

H + Y n iZ + c 2 s k 2 a- 2 n = -{{Hex)- - T^Hex), (70) 

where e\ = e/(p + p) and 

T K = -3Hc 2 s - 2H- 1 H - 2c- 1 c s . (71) 

The perturbation equation H7U|) implies the classic result |4(| , that if there is no further entropy generation then the 
curvature fluctuation 1Z will retain the constant value it reached during inflation whilst the wavelength exceeds the 
horizon size, k < aH. 

If the entropy perturbation vanishes on small as well as large scales, then we only need to solve the cosmological 
perturbation equation for 1Z. This is a reasonable approximation during supercooled inflation, and the equation is 
often used in that context |47j . 

In warm inflation the entropy perturbations are important. Nevertheless, for the case T t T = 0, it is possible to 
obtain analytic results. It is advantageous to use the inflaton equation for approximate solutions, relying on the fact 
that the metric perturbations remain small on scales smaller than the horizon, k > aH. An analytic approximation 
to the density perturbation amplitude can be obtained by matching the classical result for 84> Eq. to the thermal 
fluctuations Eq. I|37|) at the crossing time tn when k = aH, 

/7r\V2iT 5 /2rV2 T 

(72) 

This result is analogous to the result V-r = H A /(f) 2 for supercooled inflation. 
The spectral index n s is defined by 

1 dlnP K 

Us ~ 1 = ~dh[k~' (73) 

The slow-roll equations (|15I17(I enable us to express the spectral index for the amplitude Q72JI in terms of slow-roll 
parameters, 

".-l = ;(-?+t0-\p). (74) 

The first two terms agree with reference ,26]. The (3 term shows the dependence of the spectrum on the gradient of 
the reheating term. For comparison, the spectral index for standard, or supercooled inflation, is n s = 1 — 6e + 2ry |48| . 
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TABLE I: Three models with different damping terms. Tf is the temperature at the end of inflation and g, — 100. 



model 


b 


c 


A 


00 


V(0) 


To 


Tf 


I 








2 x 1CT 15 


0.3m pi 


(4.48 x 10 14 GeV0 4 


1.43 x W 13 GeV 


3.0 x 10 13 GeV 


II 


2 





6.4 x 10" 16 


0.6m pi 


(6.75 x 10 14 GeV) 4 


4.3 x 10 13 Ge\/ 


2.1 x 10 13 GeV 


III 


2 


-1 


6.4 x 1(T 16 


0.6m pi 


(6.75 x 10 14 GeV) 4 


4.3 x 10 11 GeV 


1.1 x 10 13 GeV 



V. NUMERICAL RESULTS 



We shall consider the evolution of the cosmological perturbations and the spectrum of density perturbations in 
a particular toy model which allows us to investigate a number of interesting phenomena. The toy models include 
some basic features of consistant models of warm inflation, such as [LUllfljl. but they are simplified in order to isolate 
particular effects. 

The same potential for the inflaton field will be used throughout. The thermal corrections to the potential are taken 
to be insignificant. Three different types of damping term will be considered. We hope to compare the perturbation 
spectra for different potentials and potentials with large thermal corrections in later work. 

The potential we choose is 

1/(0) = \\{tf- 4>lf. (75) 

The important features of this potential are the maximum at (j) = and the minimum at 4> — 4>o- These features can 
arise in models with spontaneous symmetry breaking and have been used before in models of supercooled inflation 

m 

For the damping term T, we will take 

r = r °(£)'Gw^ (76) 

with different choices of b and c. In the first example, b = c = 0, which allows comparison with the analytic results 
of Taylor and Berera j2?|]. The second example uses 6 = 2 and c = to isolate the effects of the 4> dependence. The 
third case, b = 2, c = — 1 is the damping term first calculated by Hosoya [Bj], and used in some consistant models of 
warm inflation, for example 



A. Homogeneous solution 

The inflaton 0, entropy s and scale factor a satisfy equations 11191) . The numerical solutions shown in figure ^ plot 
the potential and the radiation density p r for the three models listed in tabled The radiation density drops rapidly 
at first due to the expansion, but rises as the solutions approach the slow-roll regime. 

The end of inflation happens when the slow-roll parameter e = 1 + r, and coincides with the equality between 
radiation and potential energy, p r = V. It is convenient to relate the times of events during inflation to the time of 
the end of inflation, tf, and so we plot N(t), where 

N = ]n(a f /a), (77) 

along the horizontal axis. 

The radiation temperatures at the end of inflation for the three models are given in table [I] These determine the 
Hubble length at the end of inflation (when the total energy density is 2p r ), 

cHJ 1 = 1.77 x lO- 2b g~ 1/2 T^ m (78) 

where T14 = Tf / (10 14 GeV). For comparison, a comoving scale which is 500 Mpc at the present epoch would have 
length 

3.59 x lO" 2 ^ 1 m (79) 

at the end of inflation. 

In the models with b — 2, the inflation continues for larger values of 4> than the model with b = 0. A consequence 
of this is that the reheat temperature is smaller in these models. This can have an important effect on the production 
of gravitinos [5(| • 
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FIG. 1: The potential and the radiation density are plotted against N for the three models described in the text. (I) Constant 
damping, (II) Damping depends on (j), (III) Damping depends on <j> and T. In all cases r ~ 100 at N = 60. 



B. Perturbations on small scales 



Our earlier discussion of thermal fluctuations ignored the effects of fluctuations in the metric and the damping term. 
The cosmological perturbation equations imply that the metric perturbation remains small and has little effect before 
the horizon crossing time tn, when k > aH. However, for an accurate treatment, and for dealing with fluctuations in 
the damping term, we must solve the perturbed Einstein equations and perturbed equation of motion. 

Our numerical code simultaneously integrates the equations for the homogeneous background and the perturbation 
equations for a given wave number k. The integration begins with the homogeneous solution alone. Integration of 
the perturbation equations only commences at the freezeout time ip, defined by the condition k — a(rif) 1 / 2 . This is 
done in order to ensure that the noise term in the Langevin equation which we saw earlier can be neglected. When 
the integration of the perturbation equations begins, the wavelength of the perturbation is much smaller than the 
horizon size, and horizon crossing occurs at a later time tn- The integration is stopped when the radiation energy 
density reaches the same value as the vacuum energy density. 

Initial conditions have to be set for the perturbed quantities at the freezeout time when the integration of the 
perturbation equations begins. The initial inflaton fluctuations were given in Eq. 1|37[) & n d Eq. I|38l) . The initial value 
of the metric perturbation is given by the perturbed Einstein equations (|58|l and l|59|) , 

a 2 

<p= -^AirGSp, (80) 

where Sp is given by equation l|51|l . Interestingly, at the freezout time, \Sp r \ « \Sp^,\ ~ \(p + p)v\. 

The time evolution of the gauge invariant parameters describing the perturbations is shown in figure^ The entropy 
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FIG. 2: Evolution of the fluctuations for two different values of the comoving wavenumber k. The three rows represent the 
three models, constant damping (top), damping depends on <j> (middle), damping depends on (f> and T (bottom). The vertical 
lines indicates horizon crossing, k — aH. 



perturbations are on the left. In the range If <t<tf, the entropy perturbation is dominated by the radiation, 

e « 5(Ts). (81) 

The entropy perturbation has a significant peak but always approaches zero after the perturbation crosses the horizon. 
In supercooled inflation the entropy perturbation is small compared to the curvature perturbation. 

Figure |3 also shows the amplitude of the curvature fluctuations for the three damping models. The amplitudes 
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FIG. 3: The spectral index is plotted for the potential \\{<j> 2 — (f>o) 2 with (I) Constant damping, (II) Damping depends on 
4>. The physical wavenumber k/a has been normalised by the hubble length cH~ l at the end of inflation (see subsect. A), 
In (0.002 Mpc -1 c/ (a f H f )) = - ln(gi /2 Ti 4 ) - 53.6, where T14 X 10 14 GeV is the temperature at the end of inflation. 

begin relatively small and grow to approach a constant value very shortly after horizon crossing. The influence of the 
entropy perturbations can be seen in a small bump on the rising edge of the plots. The final value of the curvature 
fluctuation in model I is within 5% of the analytic approximation (|72l) . 

The third model, where the damping term depend on the temperature, differs from the other two. In this case, the 
entropy fluctuation has a direct influence on the inflaton fluctuation. This causes the curvature perturbation to follow 
the entropy perturbation more closely, and and even lead to the final value of the curvature perturbation changing 
sign. (Note that the root mean square value of the the curvature perturbation is given by the modulus in this case). 
The analytic approximation for the scalar perturbation amplitude is not applicable to this case and the numerical 
results confirm this. 



C. Perturbation spectra 

The spectral index for the first two models is plotted in figure |3] The results are in close agreement with the 
analytical formula (|74|l . When the damping depends on the temperature, the spectrum is rather more complicated 
and has the oscillatory form shown in figure 01 The location of the zeros in this spectrum depend on the parameters 
To, A and 4>q. At an extremum in the spectrum, the index changes from blue (n s > 1) to red (n s < 1) with increasing 
wavenumber k. 

The observational situation is developing rapidly. The first year of cosmic microwave background anisotropy mea- 
surements made by the WMAP satellite suggest a spectral index n s < 1 51]. However, the combination of the 
microwave data with large scale structure data |52J produce a quite different conclusion, that the spectral index runs 
from blue (n s > 1) to red (n s < 1) with increasing wavenumber k. Peiris et al [S^] give their best estimate of the 
spectral index as n s — 1.2 and dn s /dlnk — —0.077 at k = 0.002Mpc _1 . 

The possibility of a spectrum which runs from blue to red is particularly interesting because it is not commonly seen 
in inflationary models, which typically predict red spectra. Examples of inflation with blue spectra do exist however. 
The hybrid inflation models, which have two scalar fields, are examples [53ll5^ . 

The numerical results show clearly that warm inflation can produce a spectrum of density perturbations with blue 
runs to red behaviour. This behaviour is seen quite generally for potentials with spontaneous symmetry breaking 
when r has the form Eq. (|76|l with 6 = 2 and c = 0, and is a possibility when 6 = 2 and c = — 1. 

VI. CONCLUSION 

The numerical results presented here show some of the possibilities which we might expect from scalar density 
fluctuations in the warm inflationary scenario. We have focused on the effects that result from various forms of the 
damping coefficient T. Our main results are: 
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with Hosoya's damping term. The spectral index over 



• The non-stochastic evolution begins at the freezeout time tp when the fluctuations in the inflaton are given by 
equations l|3TJ) and (J2HJl. 

• The spectral index when T = T{<j)) is given by three slow-roll parameters Eq. I|74|) . In this case warm inflation 
introduces only one extra slow-roll parameter. 

• The combination of 'spontaneous symmetry breaking' potentials and T = T((f>) typically leads to n s > 1 on long 
scales and n s < 1 on short scales. 

Emphasis has been placed here on methods for solving the perturbation equations rather than constructing a 
realistic inflaton potential. However, the equations and the numerical integration have been set up in a way which 
makes them easily adaptable to different models. The inflaton potentials may include thermal corrections and different 
forms of the damping terms may also be taken into account, for example using the friction coefficients calculated in 
reference |lalln ll9 | . 

A better understanding of non-equilibrium thermodynamics will enable us to calculate the friction term in the 
inflaton equation for a wider ran ge o f conditions than is possible at present. Models of warm inflation are restricted 
by consistency requirements |l3l Il5l |26( , many of which are a result of assumptions made due to the difficulties in 
calculating far from equilibrium effects. Our results suggest that warm inflation can easily produce a spectrum of 
density fluctuations that fits the observational data and that this is a direction worth pursuing. 
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